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ABSTRACT 

We study properties of dark matter halos in a variety of models which include Dark En- 
ergy (DE). We consider both DE due to a scalar field self-interacting through Ratra-Peebles or 
SUGRA potentials, and DE with constant negative w = p/p > — 1. We find that at redshift zero 
the nonlinear power spectrum of the dark matter, and the mass function of halos, practically 
do not depend on DE state equation and are almost indistinguishable from predictions of the 
ACDM model. This is consistent with the nonlinear analysis presented in the accompanying 
paper. It is also a welcome feature because ACDM models fit a large variety of data. On the 
other hand, at high redshifts DE models show substantial differences from ACDM and substan- 
tial differences among themselves. Halo profiles differ even at z = 0. DE halos are denser than 
ACDM in their central parts because the DE halos collapse earlier. Nevertheless, differences 
between the models are not so large. For example, the density at 10 kpc of a DE ~ 10 13 M Q 
halo deviates from ACDM by not more than 50%. This, however, means that DE is not a way 
to ease the problem with cuspy dark matter profiles. Addressing another cosmological problem 
- abundance of subhalos - we find that the number of satellites of halos in various DE models 
does not change relative to the ACDM, when normalized to the same circular velocity of the 
parent halo. To summarize, the best way to find which DE model fits the observed Universe is 
to look for evolution of halo properties. For example, the abundance of galaxy groups with mass 
larger than 10 13 /i _1 M Q at z *5 2 can be used to discriminate between the models, and, thus, to 
constrain the nature of DE. 

Subject headings: methods: analytical, numerical - galaxies: clusters - cosmology: theory - dark energy 



1. Introduction 

Mounting observational evidence for Dark En- 
ergy (Perlmutter et al. 1999; Riess et al. 1998; 
Tegmark, Zaldarriaga, & Hamilton 2001; Netter- 
field et al. 2002; Pogosian, Bond, & Contaldi 
2003; Efstathiou et al. 2002; Percival et al. 2002; 
Spergel et al 2003), which probably contributes 
~ 70% of the critical density of the Universe, rises 
a number of questions regarding consequences for 
galaxy formation. Traditionally, DE is described 



by the parameter w = pj p characterizing its equa- 
tion of state. The ACDM model (to = — 1) was 
extensively studied during the last decade. Mod- 
els with a constant negative w > — 1 were much 
less studied, let alone physically motivated models 
with variable w (Mainini, Maccio, & Bonometto 
2003; Mainini et al. 2003), for which no N-body 
simulation has been performed yet. Observations 
(Spergel et al 2003; Schuecker et al. 2003) limit 
the present day value of w ^> —0.8, though the 
limit has been derived for constant-w models only. 
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In the accompanying paper (Mainini et al. 
2003) we describe procedures and give approx- 
imations for different quantities encountered in 
the linear and nonlinear analyzes of fluctuations 
in models in which DE is produced by a self- 
interacting scalar field (dynamical DE). In this 
paper we use these approximations to perform N- 
body simulations of models with dynamical DE 
and to study different properties of dark matter 
halos in such TV— body simulations. For complete- 
ness we also study models with constant w = —0.6 
and w = —0.8 

Our main interest is in the models with varying 
w. These models use physically motivated poten- 
tials of scalar field and admit tracker solutions. 
We focus on the two most popular variants of dy- 
namical DE (Wetterich 1988; Ratra & Peebles 
1988; Wetterich 1995). The first model was pro- 
posed by Ratra & Peebles (1988, RP hereafter). 
It produces rather slow evolution of w. The sec- 
ond model (Brax & Martin 1999; Brax, Martin, 
& Riazuclo 2000; Brax & Martin 2000) is based 
on simple potentials in supergravity (SUGRA). 
It results in much faster evolving w. Hence, RP 
and SUGRA potentials cover a large spectrum of 
evolving w. These potentials are written as 

A 4+a 
A4+a 

V(4>) = —, — exp(47rG</> 2 ) SUGRA. (2) 

Here A is an energy scale, currently set in the 
range 10 2 -10 10 GeV, relevant for fundamental in- 
teraction physics. The potentials depend also on 
the exponent a. The parameters A and a de- 
fine the DE density parameter Qde- However, 
we prefer to use A and FIde as independent pa- 
rameters. Figure 10 in Mainini et al. (2003) gives 
examples of w evolution for RP and SUGRA mod- 
els. The RP model considered in this paper has 
A = 10 3 GeV. At redshift z = it has w = -0.5. 
The value of w gradually changes with the red- 
shift: at z = 5 it is close to -0.4. The SUGRA 
model has w — —0.85 at z = 0, but w drastically 
changes with redshift: w « — 0.4 at z = 5. Al- 
though the w interval spanned by the RP model 
covers values significantly above -0.8 (not favored 
by observations), this case is still important both 
as a limiting reference case and to emphasize that 
models with constant w and models with variable 



w produce different results even if average values 
of w are not much different. Constant w models 
have no physical motivation and can only be justi- 
fied as toy models to explore the parameter space. 
The typical values of w observed in dynamical DE 
models, however, suggest to use w = —0.8 and 
w = —0.6 for the models with constant w. 

2. Simulations 

The Adaptive Refinement Tree code (ART; 
Kravtsov, Klypin & Khokhlov 1997) was used to 
run the simulations. The ART code starts with 
a uniform grid, which covers the whole computa- 
tional box. This grid defines the lowest (zeroth) 
level of resolution of the simulation. The stan- 
dard Particles-Mesh algorithms are used to com- 
pute density and gravitational potential on the 
zeroth-level mesh. The ART code reaches high 
force resolution by refining all high density regions 
using an automated refinement algorithm. The 
refinements are recursive: the refined regions can 
also be refined, each subsequent refinement hav- 
ing half of the previous level's cell size. This cre- 
ates a hierarchy of refinement meshes of different 
resolution, size, and geometry covering regions of 
interest. Because each individual cubic cell can 
be refined, the shape of the refinement mesh can 
be arbitrary and match effectively the geometry of 
the region of interest. 

The criterion for refinement is the local density 
of particles: if the number of particles in a mesh 
cell (as estimated by the Cloud-In-Cell method) 
exceeds the level n t h re sh, the cell is split ("refined") 
into 8 cells of the next refinement level. The re- 
finement threshold may depend on the refinement 
level. The code uses the expansion parameter a 
as the time variable. During the integration, spa- 
tial refinement is accompanied by temporal refine- 
ment. Namely, each level of refinement, I, is in- 
tegrated with its own time step Aai — Aao/2 z , 
where Aao is the global time step of the zeroth 
refinement level. This variable time stepping is 
very important for accuracy of the results. As the 
force resolution increases, more steps are needed 
to integrate the trajectories accurately. Extensive 
tests of the code and comparisons with other nu- 
merical A-body codes can be found in Kravtsov 
(1999) and Knebe et al. (2000). The code was 
modified to handle DE of different types. 
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Fig. 1. — The mass function of isolated halos in the 
ACDM and RP models. The masses are found within 
virial radii. Results of all the rest of the models are in 
between these two models. At z — there is no dif- 
ference between the mass functions. The dot-dashed 
curve shows ST prediction. At z = 2 the mass func- 
tions of w 7^ — 1 models are above that of the ACDM. 

A large number of simulations were performed. 
The simulations have different sizes of computa- 
tional box, different force and mass resolutions. 
Table 1 lists parameters of all our simulations. 
This large set of simulations allows us to study 
properties of halos ranging from dwarf satellites 
to clusters of galaxies. All simulations were exten- 
sively studied. We find that in all cases the results 
are bracketed by the ACDM and the RP models. 
Normally, differences between models are not very 
large. In order to avoid too crowded plots, in most 
of the presented plots we show only results of these 
two models. 

3. Statistics of halos: power spectrum, 
mass and velocity functions 

Figure 1 shows the mass function for isolated 
halos in the RP and the ACDM models. The sim- 
ulations have the same initial phases and the same 
value as = 0.75. Thus, the differences between 
models are only due to different w(t). Remark- 
ably, at z — the mass functions are practically 
indistinguishable: a mass function has no "mem- 



0.001 p — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — c 




12 3 4 



z 

Fig. 2. — The evolution of the number density of halos 
with virial mass larger than 10 13 /i -1 Mpc for different 
models. The bottom and top symbols are for A— body 
results in ACDM and in RP models. The curves show 
ST approximations. The dashed (dot-dashed) curve is 
for w — —0.6 (SUGRA) models. There is hardly any 
difference between models at redshifts smaller than 
z = 1. At higher z the number of halos in ACDM de- 
clines faster than for other models. 

ory" of the past evolution. In this figure we show 
only two models, but all other models show the 
same results at z = 0. The mass function is well 
fitted by the approximation provided by Sheth & 
Tormen (ST, Sheth & Tormen 1999; Sheth, Mo 
& Tormen 2001; Sheth & Tormen 2002). 

At higher redshifts the situation is quite dif- 
ferent: mass functions deviate substantially. Bot- 
tom panel in Figure 1 clearly demonstrates this: 
the number of clusters with mass large than w 
3 x 10 13 /i _1 Mq is almost ten times larger in the 
RP simulation. The differences depend on mass. 
They are larger for massive clusters and much 
smaller for less massive halos. For galaxy-size ha- 
los with mass ~ 10 12 /i _1 M© the differences are 
only ~ 20%, which will be difficult to detect ob- 
servationally. 

The dependence of halo abundance with red- 
shift is further illustrated in Figure 2, where we 
study halos with mass of a group of galaxies. 
There is almost no way to distinguish models at 
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Table 1 
Parameters of simulations 
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recent times z < 1. But at z = 2 — 3 the dif- 
ferences are quite significant. We note that ob- 
servational detection of group-size halos at high 
redshifts is difficult, but feasible. We know how 
these objects should look like - almost the same as 
nearby groups. A group at high redshift should be 
more compact than a group at z = and it should 
consist of 3-10 Milky- Way size galaxies. Galaxies 
are expected to be distorted by interactions. A 
sample of few thousands galaxies can be used to 
count the number of groups. Comparison with the 
number of groups at present moment seems to be 
the way to discriminate between different models 
of DE. 

For each halo we find the density profile and 
estimate the maximum circular velocity V c i TC = 
y/GM(< r)/r. We then construct the circular ve- 
locity function of halos - the number density of 
halos with given V^i rc . The velocity function is a 
kin of the mass function, but it probes deeper in- 
side halos. For a typical halo discussed here with a 
concentration C ~ 10, the radius of the maximum 
circular velocity is about five times smaller than 
the virial radius. Figure 3 shows the velocity func- 



0.01 F 




200 400 600 8001000 



V oirc (km/s) 

Fig. 3. — The number density of halos with the max- 
imum circular velocity larger than Kirc at z = 2. We 
show RP and ACDM models. 

tion at z = 2 for RP and ACDM models. Just as 
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Fig. 4. — Power spectrum of fluctuations of dark mat- 
ter for the ACDM (full curves) and for the RP (dashed 
curves) models at different redshifts. The dot-dashed 
curves show linear spectra for the ACDM model. At 
high redshifts fluctuations in the RP model are larger 
than for the ACDM model resulting in earlier collapse 
and in more dense halos. The differences in P(k) are 
the largest at z = 2 — 3. 

in the case of the mass function, the differences be- 
tween models are larger at high redshifts. At given 
redshift the differences are larger for massive ha- 
los. Still, the velocity function brings new results. 
Even at z = 2 the mass functions are very close 
for low mass halos with virial mass « 10 12 ft _1 M Q . 
These halos have V c i rc rs 200 km/s. The velocity 
functions at that V c - uc are visibly different: RP 
model has about 1.5 times more halos. The only 
way to explain this is to have more concentrated 
halos in RP model. In the next section we will 
explore this possibility in detail. 

Figure 4 shows the evolution of the power spec- 
trum P{k) for fluctuations in the dark matter. 
The power spectrum basically follows the same 
pattern as the mass function: relatively large 
differences at high redshift, which become much 
smaller at z — 0. At z — the deviations remain 
only on small scales (k > 2). 
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Fig. 5. — Profile of the same halo simulated in differ- 
ent models. The full curve is for the ACDM model. 
The shot (long) dashed curve is for SUGRA (RP) mod- 
els. The halo has the virial mass 5 x 10 13 /i -1 M Q . For 
different DE models the density profiles are practically 
the same in the outer R > 100/i _1 kpc region. Each 
profile is well approximated by the NFW profile. 

4. Halo structure 

We start our study of halo profiles by making 
high resolution simulations of the same halo in dif- 
ferent models. The halo was initially identified in 
a low resolution run. Short waves were added to 
the spectrum of initial perturbations and the halo 
was simulated again using « 2 x 10 5 particles. In 
the ACDM model the halo has virial mass 5 x 
10 13 h~ 1 M.Q and virial radius 730/i _1 kpc. It is ac- 
curately fitted by the NFW profile (Navarro, Frenk 
& White 2002) with the concentration C vir = 7.2. 
In the RP model the virial radius is 680ft. _1 kpc - 
visibly smaller than for the ACDM halo. The RP 
halo also have larger maximum circular velocity as 
compared with the ACDM halo. Figure 5 shows 
profiles of the halo in the ACDM, RP, and SUGRA 
models. In spite of the fact that the virial radii for 
all the models are different, the density profiles in 
the outer part of the halo R > 100ft - : kpc are prac- 
tically the same: from 100ft - : kpc to 700ft _1 kpc 
the differences are less than 10%. The halos differ 
only in the central region R < 100ft -1 kpc. The 
RP halo is clearly denser and more concentrated 
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than the ACDM halo with the SUGRA halo being 
in between. This difference can be used to discrim- 
inate between the models. Yet, it will not be easy 
because the differences are relatively small: factor 
1.5 at 10/i _1 kpc. 

The RP has a smaller virial radius because the 
virial radius in the RP model is defined at larger 
overdensity (A vir ^p — 149.8 p cr ).This is the pre- 
diction of the top-hat model of halo collapse used 
to define the virial mass (Mainini et al. 2003). 
There is nothing wrong with it, but it complicates 
the comparison of density profiles and concentra- 
tions in different DE models. For example, a halo 
with exactly the same profile will have different 
virial radii and, thus, different concentrations in 
different DE models. In order to make compari- 
son of density profiles less ambiguous, we decided 
to measure the halo concentration as the ratio of 
the radius at the overdensity of the ACDM model 
(103 times the critical density) to the characteris- 
tic ("core") radius of the NFW profile. The effect 
of using the radius at the constant overdensity in- 
stead of the virial radius is relatively small. For 
typical RP halo with virial mass ~ 10 13 /i _1 M Q the 
virial radius is ~ 15% smaller as compared with 
the constant overdensity radius. 

We also study profiles of hundreds of halos 
in simulations with lower resolution. Figure 6 
shows the dependence of halo concentration on 
the mass of halos in simulations with 80ft- _1 Mpc 
box with as — 0.75. This plot shows the same 
tendency, which we found for the high-resolution 
halo: models with dynamical DE produce more 
concentrated halos. Figure 7 shows the distri- 
bution of halo concentrations for halos in mass 
range (5 — 10) x lO 13 /i _1 M . Halos with large 
deviations from NFW fits (non-relaxed halos) are 
not used. The spread of concentrations in the 
ACDM model is about twice smaller than in Bul- 
lock et al. (2001). 

Abundance of subhalos in the ACDM model is 
a known problem (Klypin et al. 1999; Moore et 
al. 1999). It is interesting to find where dynamical 
DE models stand regarding the problem. Because 
fluctuations in dynamical DE models collapse ear- 
lier than in the ACDM model, one naively ex- 
pects that the number of subhalos is also larger. 
We study the number of subhalos in a high reso- 
lution halo. The halo is simulated in the RP and 
the ACDM models. The halo with virial mass 
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Fig. 7. — Distribution of halo concentrations for ha- 
los m mass range (5 - 10) x 1O 13 /i _1 M for different 
models. Halos with large deviations from NFW fits 
(non-relaxed halos) are not used. 
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Fig. 8. — Abundance of subhalos in a halo with virial 
mass M vir = 2.4 x lO 13 /t _1 M . When normalized 
to the circular velocity of the parent halo, the ve- 
locity function is the same for both the RP and the 
ACDM models and is well approximated by the power- 
law n(> V) oc V -2 ' 75 . Vertical bars indicate the shot- 
noise errors. 

1.3 x 10 8 /i _1 M Q . The maximum circular velocity 
of the halo in ACDM (RP) model is 522 km/s 
(594 km/s). The force resolution w l/i _1 kpc al- 
lows us to resolve dwarf DM halos with circular 
velocity larger than 30 km/s. For each (sub)halo 
we measure the density profile and estimate the 
value of the maximum circular velocity. 

The number of subhalos in the RP halo is larger 
than in the ACDM halo: Inside the radius with 
the mean overdensity 103 of the critical density 
there are 87 satellites in the RP halo and 52 satel- 
lites in the ACDM model. Thus, there are a factor 
of 1.7 more satellites in the RP halo. Nevertheless, 
this large difference can be misleading because the 
circular velocity of the RP halo is larger by factor 
1.14 and halos with larger circular velocity have 
a tendency to have more satellites (Klypin et al. 
1999). In Figure 8 we plot the number of satellites 
as the function of the ratio of the satellite velocity 
to the halo velocity. Differences between the mod- 
els are very small. It is also interesting to note 
that the velocity function of the subhalos is well 
approximated by the power-law n(> V) oc V~ 2 ' 75 . 



The slope of the power is the same as for subhalos 
of Milky Way-size halos (Klypin et al. 1999). In 
other words, it indicates that the slope does not 
depend on the mass of halo and does not depend 
on the DE equation of state. 

5. Discussion and conclusions 

Models with the dynamical DE arc in infant 
state. We do not know the nature of DE. Thus, 
a great arbitrariness exists on the choice of the 
equation of state w(t). 

At first sight it seems that the situation is hope- 
less. This paper shows that this is not true: if we 
accept that w is close to -1 at z = 0, as many ob- 
servations suggest, and that w monotonically in- 
creases with redshift, dynamical models are useful 
and can produce definite predictions for proper- 
ties of halos and for galaxies hosted by the halos. 
Furthermore, the differences between rather ex- 
treme models of DE appear to be relatively small. 
In other words, one can make detailed predictions 
for properties of dark matter halos and for their 
clustering without knowing too many details of w 
evolution. Yet, the differences between models of 
DE exist and can be used to constrain the value 
and the evolution of w. In particular, distinguish- 
ing DE models by using only the value of w at the 
present time is clearly insufficient. 

The main tendency, which we find in all DE 
models is that halos tend to collapse earlier. As 
the result, they are more concentrated and more 
dense in the inner parts. Nevertheless, differences 
are not so large. For example, the density at 
10 kpc of a ~ 1O 13 M halo in a dynamical DE 
model deviates from ACDM not more than by 
50%. This, however, means that DE is not a way 
to ease the problem with cuspy dark matter pro- 
files. Nevertheless, the differences in halo profiles 
can be exploited. Denser cluster profiles in dy- 
namical DE models can be tested by both the 
weak (Bartelmann et al 2002) and especially by 
the strong gravitational lensing. Bartelmann et al 
(1998) and Meneghetti et al (2000) argue that the 
arclet statistics favors ACDM models when com- 
pared with the open CDM models. In this respect 
dynamical DE models are between the above two 
models. This problem deserves further investiga- 
tion. 

We find that the best way to find which DE 
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model fits the observed Universe best is to look for 
evolution of halo properties. For example, com- 
parison of low- and high-z (z ^ 2) abundances of 
galaxy groups with mass larger than 1O 13 /i _1 M0 
can be used to discriminate between models. Po- 
tentially, clustering of galaxies at redshifts 2 — 3 
can also be used for this. 

In this paper we mostly pay attention to the 
group-size halos with mass ~ lO 13 /i _1 M at high 
redshifts as a probe for the DE. In the accompany- 
ing paper Mainini et al. (2003) we also argue that 
abundance of clusters at intermediate redshifts can 
be used as a test for DE models. Available cluster 
samples, unfortunately, still include too few clus- 
ters at intermediate and high redshift. 

To directly investigate the cluster mass function 
at intermediate redshift with optical instruments, 
deep optical or near infrared data are used. Ex- 
ploiting this kind of data the Red-Sequence Clus- 
ter Survey (Gladders & Yee 2000) and the Las 
Campanas Distant Cluster Survey (Nelson et al. 
2002) were compiled. Taking carefully in to ac- 
count selections effects is rather hard and these 
samples include just tenth of objects. Selection 
effects are easier to handle for clusters detected 
in X rays. The ROSAT data were used to com- 
pile a number of cluster catalogs (Ebeling et al. 
1996, 2000; de Grandi et al. 1999). The most nu- 
merous sample of flux limited clusters (REFLEX: 
Guzzo et al 1999, Schuecker et al 2003b) is based 
on the ROSAT observations. It includes 426 ob- 
jects with redshifts up to z ~ 0.3. The XMM 
Survey (Pierre 2000) will add another 800 clusters 
with redshifts up to z ~ 1. Hopefully, follow- 
up optical programs will provide redshifts for the 
clusters in the catalogs. While designed for dif- 
ferent goals, REFLEX have been already used to 
constrain many cosmological parameters such as 
a%, and, together with SNIa data, it provides im- 
portant constraints on the DE equation of state 
(Schuecker et al 2003a). 

The Suniaev-Zeldovich (SZ) effect (scattering 
of CMB photons by the hot intracluster gas) is 
even more promising for detection of high-z clus- 
ters (La Roque et al, 2003; Weller et al 2002, 
Hu 2003). The shallow all-sky survey that the 
PLANCK experiment will produce will be sup- 
plemented by narrower surveys covering a smaller 
fraction of the sky, based on interferometric de- 
vices (OCRA: Browne et al 2000; SZA: Carlstrom 



et al 2000; AMIBA: Lo et al 2000; AMI: Kneisel 
2001). 

These new cluster catalogs require more exten- 
sive and detailed theoretical modeling. Confronta- 
tion of new observational data with theoretical 
predictions will be able to discriminate between 
different DE models. 

In our analysis we also address another impor- 
tant issue: the abundance of subhalos. It is well 
known (Klypin et al. 1999; Moore et al. 1999) 
that in the ACDM model the number of predicted 
dwarf dark matter satellites significantly exceeds 
the observed number of satellite galaxies in the 
Local Group. There are different possibilities to 
explain this excess. The most attractive explana- 
tion is related with the reionization of the Universe 
resulting in heating of gas in dwarf halos, which 
prevents them from becoming galaxies (Bullock, 
Kravtsov, & Weinberg 2001,?; Somcrvillc 2002; 
Benson et al. 2002). 

We find that the number of satellites of halos, at 
z = 0, in various DE models does not change rel- 
ative to the ACDM, when normalized to the same 
circular velocity of parent halo. If the reionization 
of the Universe is the solution of the problem, then 
the DE models predict an earlier reionization of 
the Universe, because the earlier collapse of dwarf 
dark matter halos requires an earlier reionization 
to avoid too many satellites at redshift zero. The 
recent WMAP results (Kogut et al 2003, Spergel 
et al 2003), can be interpreted as giving a large 
opacity for CMB photons r ~ 0.17 ± 0.04. If true, 
this requires that the reionization occurred at a 
redshift z ri ~ 13-20, which is too large for the 
standard ACDM model (Gnedin 2000). If the 
early reionization happens in the ACDM model, 
it would predict too few satellites for the Local 
Group because too few dwarf halos collapse that 
early. Models with SUGRA DE seem to be in a 
better position to fit WMAP results and, at the 
same time, the observed number of satellites: in 
fact, in this model, halos collapse at higher red- 
shifts as compared with the ACDM model. 
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